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Abstract 

The simplest purely imaginary and piecewise constant PT-symmetric po- 
tential located inside a larger box is studied. Unless its strength exceeds a 
certain critical value, all the spectrum of its bound states remains real and 
discrete. We intèrpret such a model as an initial element of the generalized 
non-Hermitian Witten's hierarchy of solvable Hamiltonians and construct its 
first supersymmetric (SUSY) partner in closed form. 
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1 Introduction and summary 

One of the most intriguing experimental puzzles encountered in contemporary physics 
is the evident absence of SUSY partners of elementary partides in nature. In the 
context of field theory this means that SUSY, if it exists, must be spontaneously 
broken. Witten [1] proposed a schematic model which, incidentally, failed to clarify 
this breakdown but, nonetheless, survived and found a number of applications within 
the so-called SUSY quantum mechanics (SUSYQM) [2]. 

In the latter formalism one introduces the so-called superpotential W(x) and 
defines the two operators 

A=^- + W(x) À= -4~ + W(x) (1.1) 
ax ax 

with the property that the two related different (so-called 'SUSY partner') potentials 
y(±) _ e = W 2 =F W may prové both exactly solvable at the same time. An easy 
explanation of this phenomenon lies in the fact that the related Hamiltonians 

H^ = -^- 2 +V^\x)-E Q (1.2) 

become inter-related, at a convenient auxiliary energy E = E , by the factorization 
rules = ÀA and = AÀ. The spectra of # (+) and are then alike 

except possibly for the ground state. In the unbroken SUSY case, the ground state 
at vanishing energy is nondegenerate and, in the present notational set-up, it belongs 
to This means that 

A^ + \ x )=0 (1.3) 

where ipli~\x) (resp. ipl^\x)), n = 0, 1, 2, . . . , denote the wavefunctions of 
(resp. if( - )). The (double) degeneracy of ^Vn+i ( x ) > ^Pn for n — 0, 1, 2, . . . is 

implied by the intertwining relationships 

AH^^H^A H {+) A = AH ( -\ (1.4) 

In the conventional setting, the Hamiltonians (1.2) are assumed self-adjoint. 
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New horizons have been opened by the pioneering letter by Bender and Boettcher 
[3] who noticed, in a slightly different context, that the latter condition H = 
might be relaxed as redundant and replaced by its suitable weakened forms. For 
our present purposes, we shall employ their proposal and in equation (1.2) allow 
complex potentials that are merely constrained by the requirement that their real 
and imaginary parts are spatially symmetric and antisymmetric, respectively [4]. 

It is not too difficult to show that the above SUSYQM factorization scheme re- 
mains unchanged under such a non-Hermitian generalization [5, 6, 7, 8, 9]. Of course, 
the relaxation of the usual condition H = W is by far not a trivial step. Formally, 
we may put W = THT with an antilinear 'time-reversal' operator T [10]. In such a 
setting, Bender and Boettcher (loc. cit. , cf. also some older studies [11] or newer de- 
velopments [12]) merely replaced T by its product with parity V and conjectured that 
the above-mentioned and physically well-motivated weakening of Hermiticity could 
be most appropriately characterized as an antilinear 'symmetry' or "PT-symmetry', 
VTH = HVT, of all the Hamiltonians in question. Equivalently [13] one may speak 
about the P-pseudo-Hermiticity defined by the relation 

H ] = VHV~ l . (1.5) 

In this paper we intend to concentrate on implementing the resulting VT- 
symmetric SUSYQM factorization scheme in the case of the 'simplest' model which 
remains 'realistic' and 'solvable' at the same time. This means that our 'initial' 
Schròdinger equation 

ip(x) = Erj}{x) (1.6) 

(where we dropped the superscript '(+)' as temporarily redundant) will contain just 
the most trivial infmitely deep square-well form 

Í+oo x < —L 
-L<x<L (1.7) 
+oo x > L 



dx 2 



+ V r {x) + iV\x) 
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for the real part of the potential and the most elementary short-range one 



r o 



x < — l 



V\x) = < 



—l < x < 



(1.8) 



+9 



< x < l, 
x > l 



l < L, 



g>0 



for its imaginary part. As a consequence of (1.7), the wavefunctions will be defined on 
a finite interval (-L, L) with a variable length 2L, on which they satisfy the Standard 
Dirichlet boundary conditions [14, 15] 



Given the background of the result obtained in [16], we derive in section 2, an 
elegant trigonometric form of the Standard matching conditions for wavefunctions at 
the discontinuities of the potential (subsection 2.1) and discuss the practical semi- 
numerical determination of the energies with arbitrary precision (subsection 2.2). 

In section 3 we address the key concern of our present paper, viz., the inves- 
tigation of the problem in the context of SUSYQM. Here the non-Hermiticity and 
discontinuities create some specific features, which are dealt with in detail. After 
deriving the superpotential and the partner potential in subsection 3.1, we construct 
the eigenfunctions of the latter and analyze the discontinuities in subsections 3.2 and 
3.3, respectively. 

Some physical aspects of our results are fmally discussed in more detail in section 4. 

2 Trigonometric secular equation 

2.1 7^-symmetric square well inside a real one 

Let us denote the four regions — L < x < —l, — l < x < 0, < x < l, l < x < L bj 
L2, LI, Rl, R2, respectively. We shall henceforth append these symbols as subscripts 
to all quantities pertaining to such regions. The complex potential V(x), defined in 
equations (1.7) and (1.8), may therefore be rewritten as 



ip(±L) = 0. 




V L2 (x) = V L1 (x) 



ig V R1 (x) = ig V R2 (x) = 0. 



(2.1) 
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The general solution of (1.6) satisfying the conditions (1.9) can be written as 

^L2(x) — A L sm[k(L + x)\ 



ip(x) 



ipLi(x) = B L cos\i(k*x) + i^íj smh(K*x) 
^ri{ x ) — Br cosh(/tx) + i-ff smh(nx) 
k ^R2{x) = A R sm[k(L - x)] 



(2.2) 



where 



k — s + it 



E = k 2 = t 2 - s 2 



g = 2st. 



(2.3) 



Here s, t and k are real and, for the sake of definiteness, are assumed positive. A 
priori, A L , B L , Cl, A r , B r and C R are some complex constants. 

On assuming that "PT-symmetry is unbroken, we obtain the conditions 



(2.4) 



from which we get 



A R = A B* = Br = B C* r = Cr = C. 



(2.5) 



The derivative of (2.2), taking (2.5) into account, reads 

( d x ipL2{x) = kA* cos[k(L + x)] 

d x ip{x) = 



d x ipLi{ x ) = k*B* sinh(/t*:r) + i^p cosh(/t*a;) 
d x ipRi{x) = kB síïúi(kx) + i j cosh(/íx) 
. d x ïpR 2 (x) = -kAcos[k(L - x)\ 



(2.6) 



Let us now match the wavefunction and its derivative at x = and impose VT- 
symmetry in the neighbourhood of the origin: 



^i(0) = ^li(0) e 



d x ip m (0) = d x ip L1 (0)eïí 



(2.7) 



This leads to 

B,CeR. (2.8) 

It now remains to match ip and d x ip at x — ±Z. Since ip is PT-symmetric, it is 
enough to impose matching conditions at x — l: 



ipR2(l) = ipRi(l) d x %jj R2 {l) = d x ip R1 (l). 



(2.9) 



This yields 

C 

Asm[k(L-l)] = B cosh(/íZ) + i— sinh(/íZ) (2.10) 

kI 

C 

—kAcos[k(L — l)] = kBsíïúi(kI) + i— cosh(«Z). (2.11) 

i 

We conclude that the final form of ip is 

I^ L2 (x) = sm[k(L + x)} 
^fli (x) = B cosh(/ta;) + ig smh(ícx) 
^(z) = Asin[/c(L - x)] 
where the complex constant A is determined by one of the equations (2.10) and (2.11), 
while the real constants B and C have to satisfy a condition obtained by eliminating 
A between (2.10) and (2.11): 

nlB{kcos[k(L — /)] cosh(fd) + nsm[k(L — /)] sinh(fd)} 

+ iC{A;cos[A;(L - /)] sinh(/t/) + Kwa[k(L - /)] cosh(/t/)} = 0. (2.13) 

We may therefore express both constants A and C in terms of B as 

A _ B Kcsc[fc(L-f)]csch(«0 (2 u 

kcot[k(L - l)] + /ccoth(«Z) 

C = i/cZB fccot [ A: ( L ~0] cotn ( /t/ ) + /t ( 2 15) 

kcot[k(L-l)] + K,coth(K,l)' 

Since, from (2.8), the left-hand side of equation (2.15) is real, the same should be 
true for the right-hand one. The resulting condition can be written as 

k 2 cot 2 [k(L - Z)] [k coth(/t/) + k* coth(/t*Z)] 

+ k cot[k(L - Z)] [k 2 + 2kk* coth(/t/) coth(/t*Z) + k* 2 ] 

+ KK*[k COth(/C*Z) + k* coth(/t/)] = 0. (2.16) 

On expressing k 2 , k and k* in terms of s and t through equation (2.3) and using some 
elementary trigonometric identities, condition (2.16) is easily transformed into 

k sin[2A;(L - /)] [s 2 cosh(2sZ) + t 2 cos(2í/)] 
— cos [2k(L - /)][s 3 sinh(2sZ) - í 3 sin(2t/)] 
+ st 2 sinh(2sZ) - s 2 t sin(2tZ) = (2.17) 

where k = \Jt 2 — s 2 . 
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2.2 Graphical and numerical determination of the energies 

The transcendental equation (2.17) has to be complemented by the constraint (2.3), 

st = X -g. (2.18) 

The couples of roots (s n ,t n ), n = 0, 1, 2, . . . , of this pair of equations define all the 
bound-state energies E n by the elementary formula 

E n = t 2 n -s 2 n n = 0,l,2,.... (2.19) 

In practice, the (s n , t n ) vàlues may be obtained as the intersection points in the (s, t) 
plane of the curves representing the roots of the transcendental equation (2.17) with 
the hyperbola (2.18). 

Before proceeding to discuss the graphical and numerical determination of E n in 
general, it is worth reviewing three interesting limiting cases of equation (2.17). The 
first one corresponds to the limit / — > L, wherein the present square well with three 
matching points reduces to the one with a single discontinuity. Equation (2.17) then 
simply becomes 

s sinh(2sL) + t sin(2íL) = (2.20) 

which coincides with equation (9) of [14] (where g is denoted by Z and L — 1). 

The second limiting case corresponds to l — > and gives back the real square well. 
Since the constraint (2.18) then disappears, we are only left with equation (2.17) 
acquiring the simple form 

sin(2A;L) = 0. (2.21) 

Its solutions are provided by the hyperbolas t 2 — s 2 = (ff) 2 , n — 1, 2, ... , where the 
n = value is discarded because no acceptable wavefunction can be associated with 
it. We therefore arrive at the well-known quadratic spectrum E\ = (ff) 2 , n — 1, 
2, . . . , of the real square well. 

The existence of the third special limiting regime is connected with the bounded 
nature of our imaginary barrier (1.8). In the language of perturbation theory this 
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means [15] that the influence of this barrier on the vàlues of the energies (2.19) 
weakens quickly with the growth of the quantum number n. At the higher excitations, 
as a consequence, the n— dependence of the energies will not deviate too much from 
the / — > rule E n ~ n 2 3> 1. In the other words, the growth of n will imply the 
growth of t n ~ n 3> 1 and the decrease and smallness of the roots s n = g/(2t n ) <C 1. 
In this regime, we may imagine that k = t a/1 — s 2 /t 2 = t — s 2 / (2t) + 0(s 4 /t 3 ) = 
t — g 2 /(8t 3 ) + 0(l/n 7 ) so that the six components of our quantization condition (2.17), 
viz., 

s 2 k sin[2A;(L - l)] cosh(2s/) + t 2 k sin[2A;(L - l)] cos(2íZ) 
- s 3 cos[2A;(L - l)] sinh(2s/) + t 3 cos[2A;(L - /)] sin(2í/) 
+ st 2 sinh(2s/) - s 2 t sin(2í/) = 

may be characterized by their asymptotic sizes 0(l/n), 0(n 3 ), 0(l/n 4 ), 0(n 3 ), O(n ) 
and 0(l/n), respectively. Once we omit all the negligible 0(1/ n) terms and insert 
s = g/(2t) whenever necessary, we arrive at the thoroughly simplified approximate 
secular equation 

sin(2A;L) + g + O ^ = 0. (2.22) 

Its roots are easily found, 

, ... 2q 2 lL 2 m ( 1 \ , 

k = k n = —+ -1 r 1 ^- + ° H ' 2 ' 23 
2L 7r n \ n J 

and give 

i.e., a nice and elementary approximate energy formula for all the highly excited 
states. 

In the general case, the bound-state energies (2.19) of our model are determined 
from the simultaneous solutions of equations (2.17) and (2.18). Although the former 
is transcendental, one of its roots is quite obvious, namely s = t. When we realize 
that this implies k = and substitute the solution into equations (2.12) - (2.15), we 
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obtain a vanishing wavefunction. This is in accordance with an insight provided by 
the Hermitian limit g — > or l — > 0. 

The other solutions of (2.17) can be found numerically and graphically. As we can 
see in figure 1 where we work with re-scaled length units in which L — 1, they form 
semi-ovals in (s, t) plane. We can observe the absence of robustly real energy levels, 
i.e., levels remaining real for any value of g, which played their role in [16]. 

The locally decreasing character of the semi-oval màxima could cause a complexifi- 
cation of higher energy pairs while the lower pairs would remain real. In other words, 
the semi-oval màxima might be decreasing faster then the hyperbola (2.18). This 
race in decrease can be judged easily when we use a hyperbolic coordinate system. 
As shown in figure 2, in this setting, the màxima prové to increase monotonically 
while the hyperbola is represented by a horizontal straight line. Consequently, our 
model preserves a sequential merging of the energy levels. The critical value g c of 
the coupling constant g, for which the two lowest energy levels merge together, is of 
high importance. It is the boundary of exact 'PT-symmetry, which we consider to be 
physically relevant and assumed in deriving equation (2.17). For a higher value of g, 
the wavefunction PT-symmetry would be broken. 

We found g c for various vàlues of the parameter /. Since g c rises rapidly as l — > 0, 
we present its vàlues in combination of graph and table (see figure 3 and table 1). 
As the parameter l approaches zero, g c tends to infinity and the semi-oval màxima 
run to infinity as well. As explained in subsection 2.1, equation (2.17) then provides 
the bound-state energies of the real square well. On the other hand, for l — > L = 1, 
we get back the critical coupling g c ~ 4.4753, previously obtained for the square well 
in [14] and [17]. 

3 The SUSY partner potential 

The purpose of the present section is to construct and study the SUSY partner of 
the square- well Hamiltonian H^ + \ defined in equation (2.1), in the physically-relevant 
unbroken "PT-symmetry regime, corresponding to g < g c . 
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3.1 Determination of the parameters 

Identifying with the square-well potential (2.1), i.e., vj^\x) = 0, vj^\x) = —ig, 

Vri\ x ) = ' l 9i ^m\ x ) = an d Eq = kl = t'l — s'l = —Kq + ig, we obtain for the 

superpotential and the partner potential the results 

W L2 {x) = k teni[k (x + x L2 )\ 
Wli(x) = -/tgtanh[/to(x + x L i)] 
W R1 (x) = -K teaïh[K Q (x - x R1 )] 
k W R2 (x) = A; tan[A;o(a; - x R2 )\ 



W(x) = < 



(3.1) 



V<-~\x) = 



and 

' V { L2 \x) = 2k 2 sec 2 [h{x + x L2 )} 
VÍ7\x) = -2k* 2 sech 2 [4(x + x L1 )} - ig 
v m\ x ) = -2/t o sech 2 [/t (x - x R1 )\ + ig 
< V m\ x ) = 2k o sec 2 [k (x - x R2 )} 
respectively. Here xl 2 , xli, x R \ and x R2 denote four integration constants. 

We now choose xl 2 and x R2 as 

XL * = L + W XR2 = L ~2T 



(3.2) 



(3.3) 



to ensure that V^ 2 ^ and blow up at the end points x = —L and x = L. This is 
in tune with [9]. We thus get 

V^\x) = 2k 2 csc 2 [k (x + L)] V^\x) = 2k 2 csc 2 [k (x - L)]. (3.4) 

Observe that for the superpotential, W L2 (x) and W R2 (x) also blow up at these points: 



W L2 (x) = -k cot[k (x + L)] W R2 (x) = -k cot[k (x - L)}. 



(3.5) 



Let us next consider the unbroken SUSY condition (1.3), where according to (2.12) 
the ground-state wavefunction of is given by 



4m( x ) = CH = ^HMi-i)] 



,(+) 



(+) 



(3.6) 



4m( x ) = 4li(- x ) =5j +) cosh(« a;)+i-^ r siiih(«oa;). (3.7) 

KqI 

Note that the superscript '(+)' is appended to the wavefunction and the coeficients 
to signify that we are dealing with Hamiltonian H^ + \ It is straightforward to see that 



(+) 
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equation (1.3) is automatically satisfied in the regions R2 and L2 due to the choice 
made for the integration constants xr%, xl2 in equation (3.3). On the other hand, in 
the region Rl we find a condition fixing the value of xri, 

+ an U*<r \- ÍCÍ ° +) - k ° COt t fc 0( L ~ 0] COth(K l) + Kp 

tanh^J - — p - fcocotML _ z)] + ^ )Coth(Ko0 (3-8) 

where in the last step we used equation (2.15). A similar relation applies in LI, thus 
leading to the result 

x L1 = x* R1 . (3.9) 

Note that in contrast with the real integration constants Xr 2 , X L2, the constants 
xri and xli are complex. Separating both sides of equation (3.8) into a real and an 
imaginary part, we obtain the two equations 

sinhXcoshX N r 



cosh 2 X cos z Y + sinh 2 X sin Y D 
sin Y cos Y N l 



(3.10) 
(3.11) 



cosh 2 X cos 2 Y + sinh 2 X sin 2 Y D 
where we have used the decompositions /t = s + i^o> x m — x m + 1x ri, k o x ri — 
X + iY, implying that 

X = s x r R1 - t x l R1 Y = t x r R1 + s x R1 (3.12) 

and we have defined 

N r = {-s 2 cos[2A;o(L - Z)] + í 2 } sinh(2s Z) + k s sin[2k {L - Z)] cosh(2s Z) (3.13) 

N* = { s 2 - t 2 cos[2k {L - Z)]} sin(2í Z) - Mo sin[2A;o(L - Z)] cos(2í Z) (3.14) 

D = {- s 2 COS [2A; (L - /)] + í 2 } cosh(2s Z) + {4 ~ 4 cos[2k (L - Z)]} cos(2í Z) 

+ À; sin[2À;o(L - Z)][s fflnh(2s Z) + í ffln(2í Z)]. (3.15) 

Equations (3.10) and (3.11), when solved numerically, furnish the vàlues of both the 
parameters x r R1 and x R1 . 

One may also observe that the resulting superpotential W{— x) = — W*(x) and 
partner potential V^(— x) = V^*(x) are PT-antisymmetric and PT-symmetric, 
respectively. 
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= o vü 2 


tòà(-l) 













3.2 Eigenfunctions in the partner potential 

On exploiting the first intertwining relation in (1.4), the eigenfunctions ipfc~\x), n = 0, 
1, 2, of H(~ï can be obtained by acting with A on ip^^x), subject to the 
preservation of the boundary and continuity conditions 

(L) = (3.16) 
d^ítH) = d^íti'l) (3-17) 
9xtòl(0) = ^íà(0) (3- 18 ) 

Application of .4. leads to the forms 

= Cl-Uiírsm[^ +1 (L + x)] 

x {fc n+1 cot[fc n +i(L + 2;)] - k cot[k (L + x)]} (3.20) 
tó-ifa) = CÍli sinh K+i^){<+i - tanh[/t*(a; + x^)] coth« +1 a;)} 
+ ^ï^sinh(< +1 *) 

x coth« +1 x) - /t* tanh[«*(x + x* R1 )}} (3.21) 

tóàfa) = Cim sinh(/í n+1 x){/í n+ i - /t tanh[/t (x - x ra )] coth(/t n+1 a;)} 

iC (+) 

+ c Íri —^7 sinh{K n+1 x) 

K n+ iL 

x {k u+ i coth(/c n+ iz) - /t tanh[/t (a; - aju)]} (3.22) 
tó 2 (*) = C^ 2 A^ lS m[k n+1 (L-x)\ 

x {— ^n+i cot[k n+1 (L — x)] + fc cot[k (L — x)]} (3.23) 

where C^~ L2 , C^ L [, C^ ± , C^ 2 denote some complex constants and equation (3.9) 
has been used. It can be easily checked that the boundary conditions (3.16) are 
automatically satisfied by these eigenfunctions. It therefore remains to impose the 
continuity conditions (3.17) - (3.19). 

Let us first match the regions LI and Rl at x — 0. The continuity conditions 
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(3.18) yield the two relations 



C 



(-) 

nRl 



C ím { B< n+A K l+i - KoSecti 2 (K x R1 )} + :r y ±i /ío tanh(/t :r m ) 



1Ly n+l 



(-) 

nLl 



-5^ ) 1 /t*tanh(/í*a;^ i ; 



ir (+) 



(3.24) 



iC, 



(+) 



n+l 



iC 



(+) 

n+l 



«o tanh(KQ^iïi) f ■ (3.25) 



Since equations (3.8) and (2.3) provide the two constraints 



Kotanh(KoXRi) = — k$ tanh(/tQX^ 1 ) 



*2 2 
K n+1 K n+1 



*2 2 



(3.26) 
(3.27) 



equations (3.24) and (3.25) are compatible and lead to the condition 



u nRl — °nLl- 



(3.28) 



Considering next the matching between Rl and R2 at x — l, we obtain from 
equation (3.19) the two conditions 



c ÍRi{kn+iCot[k n+1 (L - l)] + K tanh[/to(/ - x R1 )]} 
= c í~À{k n +i cot[k n+1 (L - /)] - k cot[k (L - /)]} 



(3.29) 



C í~m («íl+i - K l + totanh[/t (/ - x m )\{k n+1 cot[k n+1 (L - ï)\ 

+ /t tanh[/t (/ - Xfli)]}) 
= C^ffi (fco _ ^n+i - fc o cot[fco(L - Z)]{fc„+i cot[£; n +i(L - ï)\ 

- Jfeocot[jfeo(i - Z)]}) (3.30) 

after making use of equations (2.14) and (2.15) to eliminate and C^. 

Equations (3.29) and (3.30) both yield the same result 



M-) _ r (-) 



(3.31) 
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due to the two relations 



K tanh[K (l — xri)] = —k cot[ko(L — l)] (3.32) 

and 

K n+i ~ K o = k-o ~ k-n+i (3.33) 

deriving from (3.8) and (2.3), respectively. 

Since a result similar to (3.31) applies at the interface between regions L2 and LI, 
we conclude that the partner potential eigenfunctions are given by equations (3.20) 
- (3.23) with 

Cnil = cIli = C^jji = C^jjjj = Cl \ (3.34) 
Such eigenfunctions are PT-symmetric provided we choose imaginary: 

C$T» = -C<->. (3.35) 

3.3 Discontinuities in the partner potential 

In subsection 3.1, we have constructed the SUSY partner V^\x) of a piece-wise 
potential with three discontinuities at x — — /, and /. We may now ask the following 
question: does the former have the same discontinuities as the latter or could the 
discontinuity number decrease? We plan to prové here that the second alternative 
can be ruled out. 

For such a purpose, we will examine successively under which conditions V^\x) 
could be continuous at x = l or at x = and we will show that such restrictions would 
not be compatible with some relations deriving from the unbroken-SUSY assumption 
(1.3). Observe that we do not have to study continuity at x — — l separately, since 
V^~\x) being PT-symmetric must be simultaneously continuous or discontinuous at 
x = — l and x — l. 

Let us start with the point x — l. Matching there vj^\x) and vj^\x), given in 
equations (3.2) and (3.4), respectively, leads to the relation 

-2/tg sech 2 [/t (/ - xm)] + Í0 = 2^ csc 2 [k (L - /)]. (3.36) 
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On using (3.32) and some simple trigonometric identities, such a relation can be 
transformed into fcg = — Kq + |ig, which manifestly contradicts equation (2.3). Hence 
continuity of V^~\x) at x = l is ruled out. 

Consider next the point x = 0. On equating vj^\o) with V^[\o) and employing 
(3.2) and (3.9), we obtain the condition 

— 2kq sech 2 (/t 2;m) + l 9 = — 2/íq 2 sech 2 (/tQX^ 1 ) — ig. (3.37) 

Equations (2.8) and (3.8) then yield the relation — k% + |ig = — /íq 2 — |ig, which 
contradicts equation (2.3) again. Continuity of V^~\x) at x = is therefore excluded 
too. 

We conclude that under the simplest assumption of unbroken SUSY with a factor- 
ization energy equal to the ground-state energy of H^ + \ the partner potential V^~\x) 
has the same three discontinuities at x — — /, and l as V^ + \x). 

4 Discussion 

Among all the PT-symmetric models, field-theoretical background explains the last- 
ing interest in the purely imaginary long-range model V(x) = ix 3 [18, 19] and its 
generalizations V(x) = x 2 (ix) s with the imaginary part V l (x) exhibiting, at any 

5 G [0,2), a characteristic 'strongly non-Hermitian' (SNH) long-range growth in 'co- 
ordinate' x E R. Up to the harmònic oscillator at S — 0, all of the latter SNH 
PT-symmetric models are only solvable by approximate methods. Still, rigorous 
proofs exist showing that their spectra are all real [19]. 

By rigorous means, the reality of the spectrum has also been shown for many 
other PT-symmetric potentials V. Some of them turn out to be exactly solvable 
[20, 21, 22], and those for which V l (±oo) = may be called 'weakly non-Hermitian' 
(WNH). Their WNH character is reflected not only by a less explicit influence of the 
imaginary part of the potential upon the spectrum, but also by the existence of SUSY 
partners [8, 21, 23] which, in some special cases, may be real and Hermitian [5, 21]. 

In the light of similar observation one might feel tempted to perceive WNH models 
as 'partially compatible' with our intuitive expectations. This impression may be 
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further enhanced by noticing that another exactly solvable model, viz., the typical 
WNH spiked form of the 5 = harmònic oscillator, as described in [24], proved of 
particular interest in the SUSYQM context as well [6, 23]. 

Potentials V(x) with shapes that are piece-wise constant may be considered 
equally exceptional. All of these square-well-type models with forces located inside a 
finite interval (-L, L) may be easily classified by the number of their discontinuities. 

The simplest nontrivial non-Hermitian square-well potential must have at least 
one discontinuity (= matching point at x — 0). While the real part of this V is just 
a trivial shift of the energy scale, it may be kept equal to zero. Then, the non-zero 
strength Z of the spatially antisymmetric and purely imaginary V is the only free 
(real) parameter of the whole model with SNH features [9, 14]. Its "PT-symmetry 
remains unbroken in an interval of Z G (— Z cr i t , Z cr i t ) while its ground-state energy 
becomes complex beyond Z crit ps 4.48 (in Standard units h = 2m = 1 [14, 17]). 

It is known that some of these features are genèric [15]. Quantitat i vely, their 
occurrence has also been confirmed for the twice-constant SNH model V with two 
discontinuities [16]. Qualitatively, all of these observations facilitate the applicabil- 
ity and physical interpretation of the piece-wise constant models significantly [25], 
especially because the numerical vàlues of the maximal allowed couplings prové to 
be, in general, quite large. This allows us to guarantee the (necessary) reality of the 
energies by keeping simply our choice of Z safely below this maximum. 

The family of WNH square-well models may only start at the piece-wise potential 
with three discontinuities. In our present study of such a model it was important to 
demonstrate the parallelism of its properties with the exact solutions of the smooth 
complex potentials of similar shapes [6]. 

The most obvious parallel lies in the observation that a key formal feature of the 
SUSY partners # (±) is that they may remain both non-Hermitian and PT-symmetric. 
Of course, the parity V cannot define the positive-definite norm [12, 15, 26, 27]. A 
consistent physical interpretation of the similar non-Hermitian models was recently 
agreed (cf., e.g., [28]) to lie in the existence of a new metric-like operator V(+) > 
which is positive defmite. This Hermitian operator may be assumed to play the role of 
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the 'physical' mètric [29]. Tiris means that once our equation (1.5) is satisfied by the 
old Hamiltonian and by the new, positive-definite mètric V(+), we may declare the un- 
derlying quantum Hamiltonian quasi-Hermitian, leading to the Standard probabilistic 
interpretation of the theory (cf. the recent discussions of some related subtleties in 
[30]). Against this background our attention has been concentrated upon the feasi- 
bility of bound-state construction in a model with a phenomenologically appealing 
shape of the potential. 

A couple of consequences may be expected. Our model may open the way towards 
addressing one of the most difficult problems encountered in PT-symmetric quantum 
mechanics [27], viz., the control of a possible instability of the spectrum reality [16, 31]. 
Indeed, due to the pseudo-Hermiticity property (1.5) of our Hamiltonians H, the 
energies need not be real (i.e., observable) in principle [13]. 

Our WNH model may be also characterized by the simplicity of the bound-state 
wavefunctions. This allowed us to construct the superpotential yielding access, rather 
easily, to the Witten-type SUSY hierarchy. In this regard the compact form of our 
trigono mètric secular equation was welcome and particularly important, especially for 
any future projects trying to connect the mathematical PT-symmetry with physical 
phenomenology. 

In such a perspective, the most challenging mathematical problems attached to the 
non-Hermitian models descend from the reality of their exceptional points [32] . The 
simplest solvable models of the square-well type seem to offer a transparent laboratory 
for their study since the indeterminate auxiliary pseudo-metric V coincides with the 
common parity. 

In the context of physics, the phenomenological appeal of all the piece-wise con- 
stant anàlogues of the purely imaginary cúbic force represented a strong motivation 
for the systematic constructions of the positive-definite mètric operators of [29] (cf. 
also [12, 13, 25]). In particular, the highly appealing factorized form V( + ) = CV > 
of these mètric operators has been used and, for physical reasons, the factor C itself 
has been called 'charge' (cf. [28]). For all the models with relevance in field theory 
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(like V ~ ix 3 ), the constructions of C were shown feasible by WKB and perturbative 
methods [33]. 

In comparison, the solvability of all the simpler models facilitates the construc- 
tion of C (called, usually, quasi-parity in this context [20, 24, 26, 34]). An interesting 
energy-shift interpretation of the quasi-parity (which is a new symmetry of the Hamil- 
tonian) emerged in the strongly spiked short-range model considered in [35]. 

After we return to the square-well models, the quasi-parity or charge operator C 
may be constructed in the specific form which differs sufficiently significantly from the 
unit operator just in a finite-dimensional subspace of the Hilbert space [15, 16, 25]. 
This is one of the most important mèrits of this class of models. It seems to open a 
new inspiration for a direct physical applicability of non-Hermitian models whenever 
their spectrum remains real. 
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Figure captions 



Figure 1: Solutions of (2.17) form the semi-ovals. Their in- 
tersections with the hyperbola 2st = g determine energy levels 
E = k 2 = t 2 — s 2 of the system. Here g = 650 and / = 0.04. 

Figure 2: The previous picture (Fig.l) in [ts,k] plane, where 

k = y/t 2 — s 2 . We set g = 650 and / = 0.04 again. 

Figure 3: Fifty vàlues of critical couplings g c , increasing 
rapidly as l decreases, / — > 0. 

Table captions 



Table 1: Numerical vàlues of g c in dependence on the parame- 
ter l. The table suggests that the critical coupling grows faster 
than l/l for small /. 



Table 1 



l 


1.00 


0.70 


0.50 


0.40 


0.30 


0.20 


0.10 


0.01 


0.001 


9c ~ 


4.4753 


4.8129 


6.4364 


8.6011 


13.426 


27.273 


95.832 


9895.4 


486950 
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